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Using one of the methods of quantum nonequilibrium statistical physics we have investigated the 
spin transport transverse to the normal metal/ferromagnetic insulator interface in hybrid nanos¬ 
tructures. An approximation of the effective parameters, when each of the interacting subsystems 
(electron spin, magnon, and phonon) is characterized by its own effective temperature have been 
considered. The generalized Bloch equations which describe the spin-wave current propagation in 
the dielectric have been derived. Finally, two sides of the spin transport ’’coin” have been revealed: 
the diffusive nature of the magnon motion and magnon relaxation processes, responsible for the spin 
pumping and the spin-torque effect. 


I. INTRODUCTION 


One of the central issues of spintronics is the generation 
and control of spin currents in solids. There are different 
methods to implement the spin current: optical, mag¬ 
netic, and using an electric current. The latter is espe¬ 
cially important for use in devices^ when spin-polarized 
charge carriers are injected from a ferromagnetic ma¬ 
terial to a non-magnetic. At that, in the nonmagnetic 
material at the spin diffusion length, spin accumulation 
arises. With an external perturbation usually affecting 
the kinetic degrees of freedom, the interaction between 
translational (kinetic) and spin degrees of freedom plays 
the main role in the formation of the spin response to 
the external perturbation. Combined electric dipole res¬ 
onance can serve as an example of such an effect. In this 
case, the interaction of conduction electrons and an al¬ 
ternating electric field results in resonance at the Zeeman 
frequency^. Another example of such a response is the 
spin Hall effect (SHE)&i that is exhibited as a spin cur¬ 
rent perpendicular to both the normal current and the 
spin accumulation. There also exist mechanisms of in¬ 
teraction with external fields whose energy is simultane¬ 
ously transferred to both electronic subsystems (kinetic 
and spin)£. It has been found that thermal perturba¬ 
tions may also cause the spin effects to occur. The first 
effect that has opened a new direction in spintronics - 
the influence of thermal perturbations on spin effects, 
is the spin Seebeck effect (SSE)&£. It has been turned 
out to be inherent in conducting crystals of NisiFeig. 
Afterwards, the SSE could be observed in various mate¬ 
rials such as ( Ga , Mn)As — semiconductors and metallic 
Co^MnSi ferromagnets^. Besides, later the spin Nerst 
effect (or the thermal spin Hall effect), the spin Peltier 
effect, and others have been discovered^— . As noted 
in the worbiS., there is much common between the spin 
effects implemented in an electric or an inhomogeneous 
temperature fields. So, in spintronics, studying the inter¬ 
action between charge and spin currents, the new direc¬ 
tion appears - spin caloritronic a 16 : 17 . As far back as in 
the late 20th century, a few theoretical aspects of these 
interactions were discussed^. 


Studying of the SSE in a non-conducting magnet in 
the system of a nonmagnetic conductor / magnetic insu¬ 
lator (NM/FI) La Yi Fe 5 Oiz^ has shown that this effect 
cannot be described by standard approaches as regards a 
description of thermoelectric effects^. As distinct from 
conducting crystals where the transfer of the spin angular 
momentum is due to band charge carriers, the spin See¬ 
beck effect can be realized in non-conducting magnetic 
materials through the excitation of a localized spin sys¬ 
tem. In this case, the excitations (magnons) underlying 
the spin-wave current causes the transfer of the angular 
momentum. Thus, conducting crystals and a noncon- 
ductive magnet differ from each other in type of the spin 
current, but namely - spin-wave one. It is a new type 
of the spin current. The system of localized spin can be 
deviated from its equilibrium state in various ways. In 
the experiment^, this has been achieved by passing an 
electric current through the nonmagnetic conductor Pt , 
which results in the spin accumulation and the spin cur¬ 
rent in it. The interaction of the spin-polarized electrons 
with the localized spins at the interface (NM/FI) is ac¬ 
companied by the creation (or annihilation) of magnons, 
which in turn gives rise to a perturbation of the mag¬ 
netic subsystem. Due to small spin-wave damping, the 
spin-wave current propagates at much greater distances 
than the electron spin current. This fact makes the effect 
possible to apply in practic o ^ ^ . 

An important role in the study of thermal pertur¬ 
bations is played by a lattice (phonons). Indeed, a 
inhomogeneous temperature field may cause a devia¬ 
tion of both the localized spins subsystem (m) and the 
phonon subsystem (p) from their equilibrium state. If 
the non-equilibrium state of each subsystem is charac¬ 
terized by an effective temperature T), i = m,p, the 
non-equilibrium state of the system as a whole can be 
clearly described by a set of parameters T). The dif¬ 
ference in the effective temperatures of the subsystems 
can lead to the realization of the effect - the transfer of 
the angular momentum from the magnetic subsystem to 
the lattice or to the electron spin subsystemS S^ and vice 
versa. Thermal perturbations are also responsible for the 
drag effects, whose role is essential in a range of not very 
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high temperature s 24 ! 25 . As the authors^ believe, it is 
the thermal perturbations which produce the observable 
anomaly in the temperature dependence of the SSE. 

The spin effects in non-conductive magnetic materi¬ 
als under thermal perturbations were studied in several 
paper s 22 ! 27- —. In this case, the theoretical description 
boiled largely down to the consideration of the evolu¬ 
tion of the localized moment subsystem. Given thermal 
fluctuations derived from the fluctuation-dissipative the¬ 
orem, the authors of the works mentioned above modeled 
the localized moments’ dynamics by the phenomenologi¬ 
cal Landau-Lifshitz-Gilbert equation. As to the spin den¬ 
sity dynamics in a nonmagnetic material, it is described 
by the Bloch equation with phenomenological spin re¬ 
laxation frequencies. A spin current that was injected 
into the non-magnetic material ( NM ), was calculated as 
the average value of the rate of change of the spin den¬ 
sity in a nonmagnetic material. Spin-charge kinetics in 
the systems ( F/N ), ( F/N/F) under a temperature gra¬ 
dient was thoroughly covered in the paper— in terms of 
a simple phenomenological model and the simplest ap¬ 
proximations. 

The study of the spin-thermal effects requires a correct 
description of the thermal perturbation. There are sev¬ 
eral possible scenarios accounting for such perturbations. 
For example, the reaction of weakly nonequilibrium sys¬ 
tems on a disturbance of the thermal type can be repre¬ 
sented as Fourier transform of time correlation functions 
of appropriate flux operators over the statistical equi¬ 
librium state of the system— The structure of admit¬ 
tances of such a type is analogous to expressions for the 
transport coefficients arising in the theory of equilibrium 
systems as the response to a mechanical perturbation. 
These expressions can be represented as an additional 
summand in the system’s Hamiltonian. That is why the 
response to a thermal perturbation is often found using 
indirect methods. Such cases require introducing ficti¬ 
tious external forces identical to thermal perturbations 
exerted on the system— Response to the thermal per¬ 
turbations can also be calculated based on the Onsager 
hypothesis^ about of the nature of fluctuation damping 
or by using the local-equilibrium distribution as an initial 
condition for the statistical operator of the system to be 
found— 

The response of weakly nonequilibrium systems to 
the thermal perturbations can be universally constructed 
through the method of the nonequilibrium statistical op¬ 
erator (NSO) and its modifications— The latter rep¬ 
resent the NSO as a functional of the local-equilibrium 
distribution. The method allows calculating the response 
of not only equilibrium but also strongly nonequilibrium 
systems described in terms of roughened macroscopic 
variables. Within this method, the kinetic coefficients 
are expressed via the Fourier transform of time correla¬ 
tion functions over the statistical distribution which de¬ 
scribes the unperturbed nonequilibrium process. Natu¬ 
rally, in this case, electron scattering needs to be taken 
into account. 


In the given paper, we consider the spin-thermal ef¬ 
fects in a normal conductor/magnetic insulator, using the 
NSO method to describe thermal perturbations. The ar¬ 
ticle is organized as follows. The first part formulates 
a model at hand, represents a Hamiltonian of the sys¬ 
tem and introduces basic operators and their microscopic 
equations of motion. The second part of the work cov¬ 
ers constructing both the nonequilibrium entropy oper¬ 
ator including the thermal perturbations of the system 
and the NSO. The third section of the article focuses 
on an analysis of macroscopic equations describing the 
spin-thermal effects. 

II. THE HAMILTONIAN 

Our model consists of a normal conductor (NM) and 
a ferromagnetic insulator (FI). We hold that, in the first 
conductor, the spin accumulation takes place. It can be 
obtained and implemented in various ways, for example 
by using the spin Hall effect. The conduction electrons of 
the normal metal are coupled with the localized spin sub¬ 
system of the ferromagnetic insulator through exchange 
interaction. The inelastic scattering of the electrons by 
the localized spins, accompanied by magnon emission or 
absorption, alters the electron spin orientation and un¬ 
balances the localized spin system. Consider the magnon 
scattering mechanism by the example of the magnon- 
phonon interaction. We assume the ferromagnetic in¬ 
sulator to be in a non-uniform temperature field. Let the 
system of conduction electrons in the normal conductor 
consist of two subsystems: the kinetic and spin one. In 
this case, the first is characterized by an equilibrium tem¬ 
perature T, and the second (the spin subsystem) - by a 
temperature T s . The symbol T m designates a tempera¬ 
ture of the localized spin subsystem, T p - a temperature 
of the lattice (phonons). 

The Hamiltonian of the system (NM/FI) can be rep¬ 
resented as H = Hnm + Hpi + Hl- Here 

H NM = J dx(H k (x) + H s (x)), 

^fc(x) = ^{^,<5(x-x,)}, 

3 

H s (x) = -huj s ^2 s j ^( x “ x j)- (!) 
j 

The integration is over the volume occupied by (NM). 
Sj and pj are components of the spin and angular mo¬ 
mentum of j-th electron operators, respectively. w s = 
g s poF[/h is the Zeeman precession frequency of free elec¬ 
trons in an external magnetic field directed along the axis 
z. (g a , g, o are the effective spectroscopic splitting fac¬ 
tor for electrons and the Bohr magneton, respectively); 
{A,B} = (AB + BA)/ 2. 

H F i = J dx(H m (x) + H ms (x)) (2) 
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is the Hamiltonian of of the localized spins subsystem 
. H m (x) is the energy density operator of the magnetic 
subsystem. It involves a sum of the exchange (over the 
nearest neighbors) Hss and Zeeman energy Hg- 

H S s = —J ^2 SjS j+S , H s = ^2 Sj (3) 

j 

J is the exchange integral. w m = g m poH/Ti. H ms (x) is 
the energy density operator of interaction between the s 
and m subsystems at the interface. Then 

Hms = - Jo y: f rfxs(x) S(Rj) S(x — Rj), (4) 
j 

where Jq is the exchange integral, S(Rj) being the op¬ 
erator of the localized spin with the coordinate Rj at 
the interface. The integration in (j2j) is over the volume 
occupied (FI). 

Hl is the lattice hamiltonian 

Hl = J dx(H p (x) + H pm (x)) 

H p (x) is the energy density operator of the phonon sub¬ 
system. H prn (x) is the energy density operator of inter¬ 
action between the localized spins and phonons. In the 
future, we are interested in the evolution of the magnetic 
subsystem, so the scattering of conduction electrons by 
phonons for simplicity omitted from consideration. 

THE ENTROPY OPERATOR 

To analyze the kinetics of the spin-thermal effects, we 
employ a scheme developed in the NSO method applied 
to the case of a small deviation of the system from the 
equilibrium Gibbs distribution po = exp{—5b}. The en¬ 
tropy Sq of the equilibrium system with the Hamiltonian 
H can be written as 

So = ‘ho + /3(-Hfc + H s — p n) + /3(H m + H p + H ms + H mp ), 

where y3 _1 = T is the equilibrium temperature of the 
system. 

In terms of average densities, to the non-equilibrium 
state of the system there corresponds the entropy opera¬ 
tor 

S(t) = $(t) + J dx{/3(H k (x) - p a (x,t)n a (x,t)) + 

+/3 s (x,t)H s (x,t) + /3Hp(x,t) + 

+/3 m(x, t)(H m (x, t) + H ms (x, t) + H pm (x, t)) = 

= S 0 + SS(t). (5) 

$(f) is the Massieu-Planck functional. /3 s (x, t), (3 m (x,t) 
are the local-equilibrium values of the inverse tempera¬ 
ture of the (s), (m) subsystems, respectively. p“(x, t) is 


the local-equilibrium value of the chemical potential of 
electrons with spin a =f, 

n“(x)=£«S(x-x?) 

j 

is the operator of the particle number density with the 
spin a. n(x) = to(x) + rr(x). 

The operator SS(t) = f dxSS(x,t) describes the sys¬ 
tem deviation from its equilibrium state. 

SS(t) = A J dx{5(3 s (x,t)H s (x) — /3 5p a (x,t) n a (x) + 
+<5/3 m (x,f)(R m (x) +H ms (x) + H pm (x))}, (6) 

AA = A— < A > 0 , < ... > 0 = Sp(. ..p 0 ). 

< ... >*= Sp(... p(t)), p(t) = exp {-S(t)}. 

Sp a (x,t ) = p a (x,t) - p, 5Pi(x,t) = Pi(x,t) - 13. 

( i = s, to). The quantities <5/r Q (x, f), 5/3 m (x, f) have the 
meaning of non-equilibrium additions to the chemical po¬ 
tential p and deviations of the magnetic subsystem tem¬ 
perature from the equilibrium temperature = T. 

In the linear approximation the deviation from the 
equilibrium, the NSO (or the density matrix) p(t) can 
be written as follows^ 

o 1 

p(t)=p q (t)+ J dtie €tl J dr plS(t+ ti,ti)pQ T po- (7) 

—oo 0 

Here p q (t) = exp{— S(t)} is the quasi-equilibrium statis¬ 
tical operator, S(t) being the entropy production opera¬ 
tor: 

S(t) = 6S(t) = ^- + ^(S(t),H}. 

A further algorithm for constructing the operator p(t) re¬ 
duces to finding the entropy production operator. Com¬ 
muting the operators n a (x), H s (x) with the Hamiltonian 
©, we find the operator equations of motion 

h a (x) = -VI na (x) + hf ms) (x), (8) 

H s (x) = -V I Hb (x) + H s{ms) (x). (9) 

Here I n «(x) = ^ {Pj > ^( x - x “)} is the flux density 
of particles with spin a. (x) determines the rate of 

change of the number of electrons with spin a due to the 
spin-flip electron scattering by the (m) -subsystem. 

Similarly, I Hs (x) = E,- ,<^( x - *“)} 

determines the flux density of the Zeeman energy. 
77 s ( S m)( x ) is the rate of change of the electrons s- 
subsystem local energy due to their interaction with 
magnons (to). 

A A( 7 )(x) = (ih)~ 1 [A\(x) , Hry]. 


( 10 ) 
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Revealing the spin index explicitly, we write the expres¬ 
sion for the spin current density J s (x): 

J s (x) = (l/2){n t (x) - n^x)} = — V I s z (x) + S( ms) (x), 

I a .(x) = (i?(x)-jj-(x))/2, 

S( ms) (x) = (h[ ms) (x) - nf m8) (x))/2. (11) 


From the expression m it follows that J s (x) contains 
two parts: collisionless J|(x), and collisional s/ ms )(x). 
The former is due to the flux of particles with different 
spin orientations, the latter is determined by the spin-flip 
scattering. 

Let us turn to the consideration of the magnetic sub¬ 
system. Using the method of Holstein-Primakov^, the 
hamiltonian of localized spins subsystem can be rewritten 
with spin-wave (magnon) variables (using to the opera¬ 
tors creation 6^ and annihilation 6k)- Expressed in terms 
of the magnon operators, the Zeeman and exchange in¬ 
teractions are of the form 


Hms J ^ ^ ^kt ^' I T 6q Uktl ^k'^k+qi 

k.k'.q 

(14) 

where (aka) are the creation (annihilation) operators 
for electrons with a certain spin value a ='f, The in¬ 
elastic part of the exchange interaction (H ms ) leads to 
the angular momentum transfer between electrons of the 
normal metal and magnons of the ferromagnetic insula¬ 
tor. Obviously, non-equilibrium of one of the subsystems 
(electronic or magnon) makes the magnon/spin current 
via the interface. At the same time, the total angular 
momentum conservation law provides the boundary con¬ 
dition for the current at the metal/insulator interface: 
the continuity of the spin current. 

Thus, the equation of motion for the magnetic subsys¬ 
tem can be written as 

Efm(x) VJff m (x) T ( x ) "f Em(mp)fx)- (lb) 


H s = -hujdNS + hw d Y, ^k &k- (12) 

k 


H ss - ~JNzS 2 + Y 2 JzS{l - 7k )6+ 6 k , 

k 

7k = (1/«)5Z exp{ik6}. (13) 

s 

Here N is the number of localized moments with the 
spin S', z is the number of the nearest neighbors, b £ 6k 
are the creation and annihilation operators for magnons 
with the wave vector k. For |k <5| << 1 , we have 

2(1 - 7k) ~ (I/ 2 ) . Le - w k = JST, s ( kS ) 2 - 

For body-centered and face-centered cubic lattices with 
a lattice constant equal to a, we have Wk = 2 J S{ka) 2 . In 
other words, the contribution of the exchange interaction 
for the magnon frequency has the same form as the de 
Broglie dispersion law for free particles of mass m* 


Wfc 



1 

2 JSa 2 ' 


Thus, regarding the magnon gas as free, we arrive at 


Hm = Y £ W b k6k- £ ( k ) 

k 


n 2 k 2 

2 TO* 


This expression can be interpreted as the sum of the ener¬ 
gies of the quasiparticles-ferromagnons having the quasi¬ 
momentum P with their own effective mass to* and the 
magnetic moment 37 . 

The exchange interaction (|4]) includes both the elastic 
electron scattering by the localized spins and the inelastic 
scattering process, with the former preserving the spin 
orientation and the latter happening with the electron 
spin-flip and the creation or annihilation for magnons. 
Using the creation and annihilation operators for elec¬ 
trons and magnons, we rewrite the inelastic part of the 
exchange interaction in the form 


Here In m (x) = —hu) m Is* (x) is the flux density of the 
magnon energy in the (m) -subsystem. The remaining 
summands on the right-hand side of the equation are 
responsible for the magnon-phonon scattering processes. 

Finally, the equation of motion for the lattice subsys¬ 
tem has the form: 

Hp(x) = -V Ih p {x.) + H p (rnp)(x). (16) 

Substituting the above equations of motion into the en¬ 
tropy production operator S(t) and integrate by parts 
the summands containing the divergence of fluxes. Hav¬ 
ing discarded the surface integrals, we write the entropy 
production operator in the form 

S(t) = A J cfac{-/3/ n a(x)V/z“(x,i) - /36//*(x,i)h/ ms) (x) + 

+/f/ m (x)V/3 m (x,t) + (/3 s (x,t) - /3 m (x, t))H s ( ms ) (x) + 
+ (/3 m (x,f)-/3)H m(pm) (x)}. (17) 

Revealing the spin index in the first two summands on 
the right-hand side of the formula m , we get 

I n « ( x )V/x“(x, t) = I s z (x) V /x s (x, t). (18) 

^“(x,t)n“ (ms) (x) = 

= /z s (x,t)^ (ms) (x) = —fi s (x, t)S£,(ms) (x). (19) 

Here /i s = jj) — /fo. It is obvious that the expression TO 
provides the boundary conditions for the collisional part 
of the spin current at the interface. The expressions (11811 , 
to imply that the heterogeneity in the spin accumula¬ 
tion distribution governs the collisionless part of the spin 
current. 

Before proceeding to constructing the equations for the 
averages, we show that the entropy production operator 
S ( t ) can be represented both through the spin accumula¬ 
tion /x s (x, t) and through the spin temperature /3 s (x, t). 
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For this purpose, we find a linear relationship between 
the deviations of thermodynamic coordinates and ther¬ 
modynamic forces from their equilibrium values. In this 
approximation, it is sufficient to put p(t) ~ Pq(t) and ex¬ 
pand the quasi-equilibrium operator p q (t) = exp{— S(t)} 
in powers of 5S(t) . Then we come up with the result 

S < H s (x) >*= - j cbc'{5p s (x.'H s (x')) 0 - 
-/35p a (x' ,t)(H s (x);n a {x')) 0 }, 

S < n(x) >*= - J dx , {5/3 s (x',t)(n(x);if s (x')) 0 - 
-/3<5/r Q (x / ,t)(n(x);n“(x / ))o}, (20) 


where 


5 < A >*=< A >* - < A > 0 , 


average the operator equation CEB for the spin current 
< J s (x) >= ( d/dt) < s z (x) >. Then we arrive at 


< J s (x) >= J dx' {-Dj* (x, x')/3 V 7 /x s (x' + 

+(/3m(x') -&(x'))7iw s £*, (ms) (x,x') + 

+ P /U (x', f + 1 1 ) L z (ms) (x, x')}, 

(25) 


where 


u 

L z,(m S )( X > X ') = j dt ' ^ (S(m S )( X ) S S(m S )( X '^'))> 

—oo 

0 

DJ, s z (x, x') = J dt’ e et ' V A (I A (x); 7J (x', t')). (26) 


(AB)o = J d\Sp{ApQ ABp 1 ^}. 

o 


Going over to the Fourier components of the spatial co¬ 
ordinates and taking into account that < n > t= < n >o, 
we obtain 


Ms(q^) = 



Mq),s*(-q))o 

(n(q),n(-q)) 0 


S/3 s (q, t). 


( 21 ) 


Equation m defines the relationship between the spin 
temperature T s and the spin accumulation p s . In the 
long-wavelength limit q = 0 and the steady state, we 
obtain 


p s = -hco s (1 - T/T s ). (22) 

Substituting the expression m into the equation for 
<5 < i7 s (x) > , we get 


Here A denotes time averaging, /3j(x, t + t') = /3j(x), (i = 
s,m); /x(x, t + t') = m(x). 

We have derived the generalized Bloch equation that 
describes the motion of the spin magnetization density 
- the spin diffusion and relaxation processes. The first 
summand on the right-hand side of the equation reflects 
the former process; the second summand is responsible 
for the latter that caused by the electron-magnetic im¬ 
purity interaction at the interface. In doing so, we have 
made an allowance for the effects of temporal and spatial 
dispersion of the spin diffusion tensor and spin relaxation 
frequency. 

Performing the Fourier transform over the coordinates 
and time and taking into account formulas (Oil) . (1251 ). we 
represent the right-hand side of expression as 


- {ql q x D2 X {q,u) + v zz (q,u)}P(q) 6 < s z (q,w) >, 

(n(q), s z (-q)) 0 


P(q) = 


(n(q),n(-q)) 0 ' 


(27) 


Here 


S < s z (q) >= 5/3 a (q,t) Huj s C zz ( q). 


C zz { q) = (s z (q),s z (-q)) 0 - 


(n(q),s z (-q))g 
(n(q),n(-q))o ' 


(23) 


The desired relationship between the correlation function 
C zz and the differential paramagnetic susceptibility y s — 
of the electron gas is given by 


u 

D2 X ( q.w) = ^j / q,ti)), 


(28) 

describes the diffusion of the longitudinal spin respec¬ 
tively along the D z z z z and transverse of the magnetic 
field, and 


d 

Xs = 7T77 9sPb < s z >= j3{g s pb) 2 C zz . (24) 


III. MACROSCOPIC EQUATIONS 

Let us construct the equations for the averages, which 
have the meaning of local conservation laws of the aver¬ 
age energy density for the (s) and (m) subsystems and of 
the particle number density. Inserting the entropy pro¬ 
duction operator into the expression for the NSO 0 , we 


o 

v z Z { q,faj) = / ^ie (e ““ )tl (^ ms) ( q);i( ms) (-q.il)), 

—oo 

(29) 

describes the relaxation of the longitudinal spin 
components^. In deriving the last formulas we have 
taken into consideration the spatial homogeneity of the 
unperturbed state of the system. At that, the homo¬ 
geneity exhibits in the vanishing of the correlation func¬ 
tions of the type (A(q),H(q')) only if (q' = — q). It 
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is not hard to show that the combination of the corre¬ 
lation functions for P( q) is proportional to the degree 
of polarization of the electron gas in the normal con¬ 
ductor: P(q) ~ (rA — n})/n. It should be pointed out 
that the resulting expressions for the transport coeffi¬ 
cients .D]'* (q,w),p* z (q,w) are suitable for both classical 
and quantizing magnetic fields and free from assump¬ 
tions about the nature of the spectrum, the the kind of 
statistics, etc. Similar expressions for the spin relaxation 
frequency are also inherent in the theory describing the 
spin-lattice relaxation of a nonequilibrium spin system 
with magnetic resonance saturation^. 

Averaging the operator equation (1151) . we have 

< -Hm(x) >= -V < > + < H m ( ms ) (x) > + 

T < ILm(mp) ( x ) > . (30) 

Or averaging time 

< H m (x) >= J dx' (x, x') V T /3 m (x') + 

+ L m(m S )( X > X ') (ft( x ') -/3m( x ')) + 

+ i m(mp)( X > X ') (/3m(x') -/3)}, (31) 

where 

o 

L m(m0( X ’ X ') = j dt ' ^ (^m(mI)( X ); H m{ml) (x',t')), 

—oo 

0 

= J dt' ett ' V A (I^ m (x); I] lm (x l ,t')). (32) 

— OO 

The first term on the right-hand side of the equation 
(HU is the energy density change in the magnon sub¬ 
system due to the temperature gradient, which in turn 
leads to the magnon flux. As in the case of the spin elec¬ 
tron subsystem, this term describes the magnon diffusion. 
The second and third terms of the expression are re¬ 
sponsible for the impact of the electronic (spin) and lat¬ 
tice (phonon) subsystems on the magnon energy change 
through their interaction. The role of the spin electron 
subsystem reduces to the generation and annihilation of 
the magnons by the inelastic spin-flip electron scattering 
at the normal metal/ferrodielectric interface. Accord¬ 
ing to C21) . this contribution is proportional to the spin 
accumulation /i s . The phonon subsystem proves to af¬ 
fect the magnon energy change in a twofold manner. On 
the one hand, the nragnon-phonon scattering processes 
make themselves felt in the energy relaxation behavior of 
the magnon subsystem, on the other hand, the phonon 
subsystem often acts as a ’’heating” of the heat/charge 
transfer processes by means of drag effects^. It should 
be emphasized that the booster role can also be studied 
by the NSO method. The energy transfer efficiency be¬ 
tween the phonon and magnon subsystems depends on 
temperatures of the appropriate subsystems and the de¬ 
gree of their mutual interaction. As can easily be seen 


from ED, a non-zero contribution to the spin-wave cur¬ 
rent is the consequence of the difference in the temper¬ 
atures of the above subsystems. The equation (13TT) im¬ 
plies that, depending upon its direction, the magnon flux 
produced by a temperature gradient can give rise to the 
angular momentum transfer from the magnon system to 
the electron system. Besides, the heat fluxes contribute 
to the emergence both of the thermally induced spin- 
torque effect^— and of the spin-torque effect generated 
by magnon s 55 ’ 56 . 


CONCLUSIONS 

Using one of the methods of quantum nonequilib¬ 
rium statistical physics (NSO), we investigated the 
spin transport in hybrid nanostructures: normal 

metal/ferromagnetic insulator. An approximation of the 
effective parameters, when each of the interacting sub¬ 
systems (electron spins, magnons, phonons) is character¬ 
ized by its effective temperatures was considered. We 
constructed macroscopic equations for the spin current 
caused by both the unbalanced spin subsystem and an 
inhomogeneous temperature field in the ferromagnetic 
insulator. We derived the generalized Bloch equations 
which describe the spin and spin-wave current propaga¬ 
tion in the system. At that, these allow for both the 
diffusive nature of the magnon motion and the magnon 
relaxation processes responsible for the spin pumping and 
the spin-torque effect. 


appendix 

Let us give the results coming from the expressions 
for the spin diffusion tensor (1281) . including the temporal 
and spatial dispersion and find their relationship with the 
conductivity tensor components. The correlation func¬ 
tions of the fluxes in the formula (l28l) involve the Fourier 
components of the isothermal Green’s functions 

+ OO 

G AB (t) = 6{-t) e et (A, B(t)) 0 = j ^ e™ G ba (lo) 

— OO 

(33) 

Differentiating (1331) over i, we obtain the chain of equa¬ 
tions 

— e — iujoj G = — <5(t)(13, A) + Gi, 

- C - iwo) Gx = —5(t)(B, A) - G 2 , 

. , (34) 

G(t) = G AB (t), G\(t) = 0(-t) e et (A, B(t)) 0 . 
G 2 (t)=0(-t ) e et (A, B(t)) 0 . 
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Here A = (ih) 1 [A, Hy], Fly is the scatterer Hamilto¬ 
nian. The formal solution of the chain is 


G A b{u>) = [Mab(w) + e - i(w - w 0 )] l (A,B), (35) 

where M(u>) = G\ G -1 is the mass operator for the 
Green’s function. 

Let us consider the frequency dispersion by the ex¬ 
ample of the longitudinal spin diffusion coefficient. Re¬ 
stricting ourselves to the Born approximation of the in¬ 
teraction with the lattice in the expression for the mass 
operator, we have 


D 


ZZ 

ZZ 


1 o 

C zz izff (w) - iu' 


(36) 


Calculating the correlation functions (A,B) o, we come 
up with 


(I£*>I%*)o = n/2m, 


_ n F_ 1 / 2 (h/T) 
~ 8 F 1/2 (^/T) ’ 


where F m (x) are the Fermi integrals. The expression 
v z z z z coincides exactly with the formula for the relaxation 
frequency v v of the electron momentum^. Thus, the 
longitudinal spin diffusion coefficient is given by 


DU = Do 


— VjJ 


2T F_ 1/2 (jji/T) 
v P F 1/2 (h/T) ' 


(37) 


The components of the spin diffusion tensor can be ex¬ 
pressed in terms of the components of the conductivity 
tensor <7jfc(q, ui), which in our notation is 



o 

J dte^-^ (P N (q),I k N (-Q,t)). 


— OO 


(38) 


Simple calculations show that 

The given work has been done as the part of the state 
task on the theme ’’Spin” 01201463330 (project 12-T-2- 
1011) with the support of the Ministry of Education of 
the Russian Federation (Grant 14.Z50.31.0025) 
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